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Abstract. We report on the results of deep narrow-band Lya and broad-band U and I imaging of the fields of two Gamma-Ray 
bursts at redshift z=2.04 (GRB 000301C and GRB 000926). We find that the host galaxy of GRB 000926 is an extended (more 
than 2 arcsec), strong Lya emitter with a rest-frame equivalent width of 71 j^jjj A. The galaxy consists of two main components 
and several fainter knots. GRB 000926 occurred in the western component, whereas most of the Lya luminosity (about 65%) 
originates in the eastern component. Using archival HST images of the host galaxy we measure the spectral slopes (f\ oc A* 9 ) 
of the two components to /3 = — 2.4±0.3 (east) and — 1.4±0.2 (west). This implies that both components contain at most small 
amounts of dust, consistent with the observed strong Lya emission. The western component has a slightly redder V— I colour 
than the eastern component, suggesting the presence of at least some dust. We do not detect the host galaxy of GRB 000301C 
in neither Lya emission nor in U and I broad-band images. The strongest limit comes from combining the narrow and U- 
band imaging where we infer a limit of U(AB)>27.7 (2a limit per arcsec 2 ). The upper limits on the Lya flux implies a Lya 
equivalent width upper limit of ~ 150 A. We find eleven and eight other galaxies with excess emission in the narrow filter in the 
fields of GRB 000301C and GRB 000926 respectively. These galaxies are candidate Lya emitting galaxies in the environment 
of the host galaxies. Based on these detections we conclude that GRB 000926 occurred in one of the strongest centres of star 
formation within several Mpc, whereas GRB 000301C occurred in an intrinsically very faint galaxy far from being the strongest 
centre of star formation in its galactic environment. Under the hypothesis that GRBs trace star formation, the wide range of 
GRB host galaxy luminosities implies a very steep faint end slope of the high redshift galaxy luminosity function. 
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1. Introduction 

In addition to the very important question of what may be the 
central engine of Gamma-Ray Bursts (GRBs, see van Paradijs 
et al. 2000 for a recent review) these enigmatic events consti- 
tute a powerful new tool for the study of the intermediate and 
high redshift galaxy population. This is mainly because GRB 
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selection of galaxies is independent of other selection meth- 
ods such as continuum flux (e.g. Steidel & Hamilton 1992; 
Steidel et al. 1996; Adelberger & Steidel 2000; Fontana et al. 
2000), gas absorption (Wolfe et al. 1986; M0ller & Warren 
1993, 1998; Leibundgut & Robertson 1998; Fynbo et al. 1999; 
Kulkarni et al. 2000, 2001; Ellison et al. 2001; Warren et al. 
2001; M0ller et al. 2001) or Ly-a emission (see below). By 
comparing the results of these very different methods by which 
to select high redshift galaxies we can better understand the se- 
lection biases inherent to each method and hence obtain a more 
complete understanding of the underlying population of high 
redshift galaxies (see also M0ller et al. 2001). 

Since the first detection of an optical afterglow of a GRB in 
1997 (van Paradijs et al. 1997) the detection of about 25 GRB 
host galaxies have been reported (Hogg & Fruchter 1999 and 
references therein; Bloom et al. 1999; Holland & Hjorth 1999; 
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Hjorth et al. 2000; Klose et al. 2000; Fruchter 2001; Vreeswijk 
2001; Castro-Tirado et al. 2001; Holland et al. 2001; Fynbo et 
al. 2001a, Fruchter et al. 2001a; Kaplan et al. 2001; Fruchter 
et al. 2001b; Bloom et al. 2001). Only for GRB 990308 has 
an unambiguous detection of the host galaxy not yet been 
reported. GRB host galaxies have been detected in a very 
broad range of redshifts from z=0.43 (GRB 990712) to z=4.50 
(GRB 000131) excluding GRB 980425 at z=0.0085 which may 
be a burst of different nature. Furthermore, the measured opti- 
cal magnitudes of GRB hosts also span a very wide range from 
V=22.8 (GRB 980703) to V«30 (GRB 980326, GRB 990510, 
GRB 000301C). This span is both due to a wide range of dis- 
tance moduli and a large spread in absolute magnitudes. In fact, 
GRB-selection of galaxies is currently the only technique that 
allows the detection of galaxies over such a wide range of red- 
shifts and absolute luminositites. Most host galaxies detected 
so far are estimated to be fainter than L* at similar redshifts 
(L* is here defined from the fit to a Schechter function to the 
observed luminosity function at the given redshift). This fact 
may be due to intrinsic extinction in the GRB host galaxies 
(Sokolov et al. 2001), but it is not known if GRB host galax- 
ies have more than average intrinsic extinction compared to 
other galaxies at similar redshifts. In general GRB host galax- 
ies are actively star-forming, i.e. characterised by restframe UV 
emission and large emission-line equivalent widths (Bloom et 
al. 2002; Fruchter et al. in prep.). Searches for sub-mm emis- 
sion from GRBs (Smith et al. 1999, 2001) show that GRB host 
galaxies are in general not similar to the obscured star bursts 
detected in deep SCUBA and ISO surveys (e.g. Ivison et al. 
2000; Franceschini et al. 2001). 

In this paper we present a study of the host galaxies 
and environments of the two z«2 bursts GRB 000301C and 
GRB 000926. The goal is to place the host galaxies and their 
galactic environments in the context of other z«2 galaxies. 
The redshifts of the two bursts have both been determined 
via absorption lines in the spectra of the optical afterglow. 
GRB 000301C was at z=2.040 (Smette et al. 2001; Castro 
et al. 2000; Jensen et al. 2001) and GRB 000926 at z=2.038 
(Fynbo et al. 2001c; Castro et al. 2001; M0ller et al. in prep.). 
The host galaxy of GRB 00030 1C is very faint with esti- 
mates of its magnitude between R=28.0 (Bloom et al. 2002) 
and R=29.7 (Levan, private communication; Fruchter et al. in 
prep.), whereas in the case of GRB 000926 Fynbo et al. (2001a) 
and Price et al. (2001) report the detection of an extended R=24 
galaxy consisting of several compact knots near the position of 
the optical afterglow. Here we use narrow-band Lya observa- 
tions to probe these two host galaxies and to search for other 
galaxies in their environments. It is presently unknown whether 
the galactic environments of GRB host galaxies are overdense 
or representative of the general field. Narrow-band Lya obser- 
vations are ideal to address this question as it has been shown to 
be an efficient way to probe the faint end of the galaxy popula- 
tion at redshifts 2<z<5 (M0ller & Warren 1993; Francis et al. 
1995; Pascarelle et al. 1996; Thommes et al. 1997; Pascarelle 
et al. 1998; Cowie & Hu 1998; Hu et al. 1998; Fynbo et al. 
1999, 2000; Kudritzki et al. 2000; Kurket al. 2000; Pentericci 
et al. 2000, Steidel et al. 2000; Roche et al. 2000; Rhoads et 
al. 2000; Malhotra & Rhoads 2002; Fynbo et al. 2001b). Lya 



selection has the advantage that it is not continuum flux limited 
and hence allows the detection of high redshift proto-galaxies, 
or Lya Emitting Galaxy-building Objects (LEGOs, M0ller & 
Fynbo 2001), that are intrinsically much fainter than those se- 
lected in continuum flux limited surveys. However, since not 
all galaxies are Lya emitters (about 30% of the R<25 z«3 
Lyman-Break galaxies are (Steidel et al. 2000)), galaxy sam- 
ples selected from Lya emission will not be complete to the 
continuum flux limit of the sample. 

The paper is organised in the following way. In Sect. 2 we 
describe the observations and the data reduction, in Sect. 3 we 
present our results and in Sect. ^ we discuss the implications 
of our results. We have assumed a cosmology where Ho = 65 
km s _1 Mpc -1 , fi m = 0.3, and J7a = 0.7. For this cosmology 
a redshift of 2.04 corresponds to a luminosity distance di um 
= 17.16 Gpc and a distance modulus of 46.17. One arcsecond 
corresponds to 9.00 proper kpc and the look-back time is 11.1 
Gyr. 

2. Observations and Data reduction 

The observations were carried out during five dark, photomet- 
ric/clear nights in May 2001 and four dark nights in August 
2001 at the 2.56-m Nordic Optical Telescope (NOT) using the 
Andalucfa Faint Object Spectrograph and Camera (ALFOSC). 
The ALFOSC detector is a 2048 2 pixels thinned Loral CCD 
with a pixel scale of 07 189. 

The fields of GRB 00030 1C and GRB 000926 were im- 
aged in three filters: the standard U and I filters and a spe- 
cial narrow-band filter manufactured by Custom Scientific Inc. 
The narrow-band filter (CS3701/45) is tuned to Lya at z=2.04 
and has a width of 45 A (corresponding to a redshift width of 
Az = 0.037 for Lya or a Hubble flow depth of 3600 km s" 1 ). 
The filter transmission curves of the CS3701/45, U and I filters 
are shown in Fig. [j]. 

In order to reduce the contribution to the total noise 
from read-out noise, the detector was binned 2 by 2 dur- 
ing the narrow-band observations, and the integration times 
for the individual exposures was set to 4000 s. The total in- 
tegration times for the GRB 000301C field was 12.2 hours 
(CS3701/45), 7.3 hours (U-band), and 3.0 hours (I-band). For 
the field of GRB 000926 total integration times were 12.2 hours 
(CS3701/45), 7.3 hours (U-band), and 3.1 hours (I-band). The 
journal of observations is presented in Table [j] The individual 
exposures were bias-subtracted and flat-field corrected using 
standard techniques. The individual reduced images were com- 
bined using a cr-clipping routine optimised for faint sources 
(for details see M0ller & Warren, 1993). The full-width-at-half- 
maximum (fwhm) of point sources in the combined images are 
0"83 (I-band), l'/25 (U-band) and 1'.'24 (CS3701/45). 

In order to calibrate the narrow-band observations the 
four HST spectrophotometric standard stars feige66, feige67, 
BD+25°4655, and BD+33°2642 (Colina & Bohlin 1994) 
were observed. We obtained well sampled Spectral Energy 
Distributions (SEDs) of these four stars from the archive of the 
HST Calibration Data Base System[] We then calculated the 

1 ftp://ftp.stsci.edu/cdbs/cdbsl/calspec/ 
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Table 1. The log of observations at the NOT, May and August 
2001. 
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Fig. 1. Upper panel: The transmission curves of the three fil- 
ters used in this study, namely the narrow CS3701/45 filter and 
the broad-band U and I filters. Lower panel: The transmission 
curve of the CS3701/45 filter. The vertical lines mark the wave- 
lengths of Lya redshifted to GRB 00030 1C (dashed line) and 
GRB 000926 (full-drawn line). 



AB magnitudes of the stars in the narrow filter as the weighted 
mean of the AB magnitudes within the filter, where the filter 
transmission was used as weight. We find AB magnitudes of 
9.96, 1 1 .08, 9.06, and 10.64 for feige66, feige67, BD+25°4655, 
and BD+33°2642 respectively. Instrumental magnitudes were 
measured in a 1 1 arcsec circular aperture. To the difference be- 
tween standard AB magnitudes and instrumental magnitudes 
as a function of airmass we finally fitted a zero-point and an 
extinction term. For counts given as electrons per second we 
find a zero-point of 21.57 and an extinction term of 0.42 mag- 
nitudes per unit airmass. The broad-band images were cali- 
brated using secondary standards from Henden et al. (2000) 
and brought onto the AB-system using the transformations 
given in Fukugita et al. (1995): U(AB) = U+0.69; I(AB) = 
1+0.43. 



3. Results 

Contour plots of the combined narrow-band images of the 
340x340 arcsec 2 fields surrounding the two GRB positions 
are shown in Fig. ||. The position of 19 candidate LEGOs (see 
Sect. 3.3 below) are indicated with circles, and the position of 7 
galaxies with a deficit of flux in the narrow filter (see Sect. 3.5) 
are shown with squares. The position of GRB 000301C is indi- 
cated with an 'x'. 



3.1. Object Detection 

In order to create a catalogue of reliable candidate emission- 
line galaxies we need robust detection and precise colours. We 
use the software SExtractor (Bertin & Arnouts, 1998) for both 
detection and photometry. SExtractor uses a detection image to 
build the isophote apertures and then measures the fluxes of ob- 
jects in the narrow, U, and I images in these apertures. Because 
the narrow band is inside the U-band, the optimal detection 
of emission line objects is obtained on a detection image cre- 
ated as a weighted sum of the combined narrow-band and U- 
band images. We found that the weights which optimized the 
Signal-to-Noise (S/N) ratio of emission line objects were 0.7 
and 0.3 for the narrow-band and U-band images respectively. 
This combined image was then convolved with a Gaussian fil- 
ter with a width similar to that of the combined seeing. For the 
SExtractor input detection parameters we used a minimum area 
of 5 (binned) pixels 1.1 a above the noise in the background. In 
our final catalogue we only include objects detected at a total 
S/N ratio larger than 5 in an isophotal aperture in the narrow- 
band. In total we found 648 objects in the two fields. Objects 
detected at the 5a levels in the narrow-band have AB isophotal 
magnitudes of approximately 25.5^ This corresponds to an ob- 
served emission line flux of 2x 1 0~ 17 erg s -1 cm~ 2 (for details 
on the photometry see Sect. |3.2| ). In a 3'.'0 circular aperture the 
5<7 limit is n(AB)=25.0 corresponding to 3.5xl0~ 17 erg s~ x 
cm -2 . The combined U-band image is 1 magnitude deeper for 
continuum sources, i.e. U(AB)=26.5 at 5<r significance in the 
isophotal aperture. 



In the following we shall denote with "n" the magnitude in the 
narrow filter. 



4 J.P.U. Fynbo et al.: Deep Lya imaging of two z=2.04 GRB host galaxy fields 




Fig. 2. The combined images of the field of GRB 000301C (left) and GRB 000926 {right) each based on 12.2 hours of narrow- 
band imaging. The positions of 19 candidate emission-line galaxies in the two fields have been marked with circles. The position 
of the host galaxy of GRB 00030 1C, for which no emission is detected, is marked with "x". The positions of 7 galaxies with a 
deficit of flux in the narrow-band filter are shown with squares. The corners of the images are lost due to vignetting. The effective 
area of each field is 27.6 arcmin 2 . East is left and north is up. The distance between the tick-marks is 1.14 arcmin. 



3.2. Photometry 



3.3. Selection of Candidate LEGOs 



To derive colours for all detected objects we measure the 
isophotal magnitudes. The isophotal apertures are defined on 
the detection image, and the same (small) aperture is then used 
for all three frames narrow, U and I. To the extent that the see- 
ing of the images in the different passbands are equal, which 
is the case for the U and narrow-band images, and if the mor- 
phologies of the objects are independent of colour the isophotal 
magnitudes therefore provide photometry appropriate for ac- 
curate colour determination. The systematic error on the U-I 
colours caused by the slightly better seeing in the combined I 
image is much smaller than the photometric errors and hence 
ignored. To estimate the total magnitudes of the objects the 
isophotal aperture is too small. Instead we use aperture pho- 
tometry with 3'.'0 and 4'.'0 apertures and determine aperture cor- 
rections based on the radial profiles of the candidate emission 



galaxies as described in Sect. 3.4 



We calculate two-sided error-bars on magnitudes as 

a(m) + = mag(flux — o-(fmx)) — mag(flux) 

a(m)~ — mag(flux) — mag(flux + cr(flux)), 

where a(flux) is derived by SExtractor. We calculate error- 
bars on colours using the maximum likelihood analysis de- 
scribed in appendix A. 



For the final selection of candidate LEGOs we use the "narrow 
minus on-band-broad" versus "narrow minus off-band-broad" 
colour/colour plot technique (M0ller & Warren 1993; Fynbo 
et al. 1999, 2000b). In order to constrain where objects with 
no special spectral features in the narrow filter fall in the di- 
agram, we calculate colours based on synthetic galaxy SEDs 
taken from the Bruzual & Chariot 1995 models. We have used 
models with ages ranging from a few Myr to 15 Gyr and with 
redshifts from to 1.5 (open squares) and models with ages 
ranging from a few Myr to 1 Gyr with redshifts from 1.5 to 
3.0 (open triangles). Higher redshifts are not relevant for our 
3700 A selected survey. For the colours of high redshift galax- 
ies we include the effect of intergalactic Lyman line absorption 
(M0ller & Jakobsen 1990; Madau 1995). In Fig. | we show the 
n(AB)-U(AB) versus n(AB)-I(AB) colour diagram for the 
calculated galaxy colours (upper panel) and for the observed 
sources in the two GRB fields (middle and lower panels). The 
open squares in the upper panel shows the colours of the model 
SEDs. The dashed line indicates where objects with broad- 
band colour equal to that of the host galaxy of GRB 000926, 
and with either absorption or emission in the narrow filter, will 
fall. Emission line objects will fall in the lower left corner of 
the diagram (due to excess emission in the narrow-filter). In 
the middle panel we show the colour-colour diagram for all ob- 
jects detected in the two fields. Several objects are found to lie 
significantly away from the locus of continuum objects. 
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Fig. 3. Calculated and observed colour-colour diagrams. Upper 
panel: Calculated colour-colour diagram based on Bruzual and 
Chariot galaxy SEDs. The filled squares are 0<z<1.5 galax- 
ies with ages from a few to 15 Gyr and the open triangles are 
1.5<z<3.0 galaxies with ages from a few Myr to lGyr. The 
dotted box contains all these calculated galaxy colours. The 
dashed line is the locus of objects having the same broad-band 
colours as GRB 000926 host galaxy and various amounts of 
absorption (upper part) or emission (lower part) in the narrow- 
band filter. Middle panel: Colour-colour diagrams for all ob- 
jects in the two GRB fields. The squares with error-bars indi- 
cate objects detected at S/N>5 in the narrow-band image. As 
expected, most objects have colours consistent with being in 
the dotted box, however, a number of objects, including the 
GRB 000926 host, are seen in the lower left part of the dia- 
gram. Lower panel: The colours of the 19 candidate LEGOs. 
The rank 1 and 2 candidates are shown with filled squares and 
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The probability that a Lya emitter will separate itself sig- 
nificantly from the pure continuum objects is a (complex) func- 
tion of its continuum magnitude, colour, and Lya emission line 
equivalent width. It is therefore necessary to define a strat- 
egy to most effciently identify the Lya emitters. One may 
decide to err on the low probability side (providing a large 
list of candidates but with fairly low confirmation fraction, 
e.g. ~30% in the work by Arnaboldi et al. 2001 and refer- 
ences therein) or on the high probability side (shorter can- 
didate list, but higher confirmation efficiency, e.g. ^100% in 
Warren & M0ller (1996), Kudritzki et al. (2000), and Fynbo et 
al. (2001b)). Both strategies have their merits and their draw- 
backs; the conservative candidate list will likely miss a num- 
ber of objects while the non-conservative candidate list in the 
end may provide a few more confirmed LEGOs, but at the 
price of a much lower efficiency during spectroscopic confir- 
mation. However, when spectroscopic follow-up is done via 
MOS (Multi-Object Spectroscopy) one can to some extent have 
the best of both strategies as the less secure candidates can be 
used for slits not covered by the secure candidates. We there- 
fore decided to make a ranked candidate list in the following 
way: We first select all objects with S/N>5 in the narrow-band 
image and an isophotal colour n(AB)— U(AB)<— 0.7. We find 
32 such objects. We subsequently rank these candidates in four 
groups, \:certain, 2:good, 3:possible, and 4:rejected. For the 
rank 1,2, and 3 groups we expect rough confirmation efficien- 
cies of ~400%, ~75% and ^50% respectively. We base the 
rank on visual inspection (to exclude spurious detections or ob- 
jects with wrong colours due to bright neighbours) and on the 
aperture photometry (the aperture colour should also be consis- 
tent with excess emission in the narrow filter). In the following 
we will only discuss the 12 and 7 candidates that fall in the 
certain and good categories (our conservative candidate list). 
The possible candidates we plan to use as fill-up objects during 
the follow-up spectroscopic observations. The colours of the 
rank 1 and 2 candidates are shown in the lower panel of Fig. [| 
We find eleven rank 1 and 2 candidate LEGOs in the field of 
GRB 00030 1C and eight in the field of GRB 000926 (including 
the host galaxy). In Table || and Table || we list the photometric 
properties of the rank 1 and 2 candidates. In Fig. || we show 
regions of size 10x10 arcsec 2 around each candidate from the 
combined narrow-band (top row), U-band (centre row), and I- 
band (bottom row) images. In Fig. || we show the distribution of 
the candidate LEGOs in each of the two fields projected on the 
sky. In the field of GRB 00030 1C (the left plot) the candidates 
are all located in the southern part of the field. In the field of 
GRB 000926 the eight galaxies are more uniformly distributed 
over the field of view. The skewed distribution of LEGOs in the 
field of GRB 00030 1C suggests a large scale structure, possi- 
bly a filament similar to the one found in the field of the z=3.04 
QSO Q1205-30 (M0ller & Fynbo 2001). Spectroscopic follow- 
up is needed to place the LEGOs in 3D redshift space. 



Below we describe how we determine Lya fluxes, Lya 
equivalent widths and star formation rates for the candidate 
LEGOs. 
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Ly a fluxes: We determine Lya fluxes, fh ya , as 

fhya = fx ■ J T , CS3701/45(A)dA/T™ s a 3 X 701/45 , 

where Tcs37oi/45 (A) is me transmission curve of the narrow 
filter, 2^^701/45 * s t ^ le to P transmission of the filter, the term 
/ T cs37oi/45(A)dA/T™ s a 3 x 701/45 is the "effective" width of the 

filter (equal to 45 A), and fx is the average specific flux of 
the source determined from the narrow-band AB magnitude, 
n(AB), as 

f _ t C _ i n -0.4(n(AB)+48.60) c 
,/A — Ju-jp — J.U -^2 

Lya equivalent widths: To estimate equivalent widths we derive 
the relation between the n(AB)-U(AB) colour and the observed 
equivalent width (W om (o6s)). The n(AB)-U(AB) colour for an 
object with specific flux /„ is given by the expression: 

_„ , / fv M?CS3701/45 {v)dv/ $ T C S3701/45 0^ 

' ° g Jf u (v)T v (v)dv/jT v (v)dv 

where Tcs37oi/45 an d Tu are the transmission curves of the 
narrow and broad-band U filters respectively. We assume that 
the specific flux f\ is a constant, fo, except for an emission line 
with Gaussian profile at Ao = 3701 A with equivalent width 
W cm (obs). Under these assumptions we have 

fx = fo + foG(X - A ) x W em (obs), 

where G(A — Ao) is a Gaussian function with area equal 1. 

Finally, we use /„ oc f\ A 2 and insert in the expression for 
n(AB)-U(AB) above to obtain the relation between W em (obs) 
and n(AB)-U(AB). Assuming instead that f v is constant leads 
to practically the same relation. 

Star formation rates: Kennicutt (1983) has derived the follow- 
ing scaling relation between the Star Formation Rate (SFR) and 
the Ha luminosity L(Ha): 

SFR = L(Ha)/1.12 x 10 41 ergs _1 

Assuming negligible dust extinction we can use the case B re- 
combination (relevant for gas that is optically thick in the H I 
resonance lines) ratio L(Lya)/L(Ha)=8.7 (Brocklehurst 1971) 
to derive L(Ha), where L(Lya) = /Lya^Trd? , and hence 
the SFR. The effects of dust will be to reduce the Lya flux 
for a given SFR (Chariot & Fall 1993; Valls-Gabaud 1993; 
Deharveng et al. 1995; see also Neufeld 1991). Therefore SFRs 
derived in this way must be considered lower limits. 

The Lya fluxes, Lya equivalent widths and SFRs for the 
candidate LEGOs appear in Table || and Table [|. 

3.4. Extendedness of Lya emission 

Previous studies have found evidence that the Lya emission is 
more extended than the continuum emission (M0ller & Warren 
1998; Fynbo et al. 2001). To test whether this is also the 
case for the Lya emitters in this study we proceed as follows. 
Because the individual sources are too faint for precise profile 
fits, we stack the isolated rank 1 candidates and measure the 




III 

: 1^ 


i i 1 


I l U 





R (arcsec) 

Fig. 4. Top left: PSF (calculated as azimuthal average) of the 
combined images of I, U and narrow. Bottom left: Narrow -band 
image azimuthal average of stacked point sources (solid line) 
and of stacked Lya emitters (LEGOs, shown with x). The 
dotted lines indicate the ±ler range. The LEGOs are clearly 
extended. Bottom right: Broad-band (sum of I and U) image 
azimuthal average of stacked point sources (solid line) and of 
stacked LEGOs (x). The dotted lines indicate the ±ler range. 
The S/N of the broad band profile is lower than for the narrow- 
band profile, but there is marginal evidence that the continuum 
emission of the LEGOs are less extended than the Lya emis- 
sion. 



radial profile on the stacked image in narrow-band, U, and I. 
We compare the profiles with the profile of the same number 
of point sources stacked in the same way. The results of this 
procedure are shown in the three panels of Fig. || The Lya 
objects are clearly extended. Furthermore, although the S/N is 
lower than for the narrow-band profile, there is marginal evi- 
dence that the continuum sources of the Lya objects are less 
extended than the Lya emission itself. 

We finally use the profiles of the point-sources and the can- 
didate Lya emitters to determine the aperture corrections for 
our 3'.'0 and 4'.'0 circular apertures relative to a 1 1 circular aper- 
ture (as used for the standard stars, Sect. ||). We find aperture 
corrections for the 3'.'0(4'.'0) apertures of 0.36(0.21), 0.27(0.15) 
and 0.23(0.14) mag for narrow, U and I respectively. 

3.5. Narrow-band absorption objects 

Roughly half of the Lyman-Break selected Galaxies have Lya 
in absorption (e.g. Shapley et al. 2001, their Fig. 7). In order 
to detect z=2.04 galaxies with absorption in the narrow-band 
filter, we run SExtractor using the combined U-band image as 
detection image, hence optimising the detection of objects that 
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are bright in the U-band. We detect 515 and 578 objects above 
S/N=5 in isophotal apertures in the fields of GRB 00030 1C and 
GRB 000926 respectively. Since the observations were tuned 
to select objects with excess emission in the narrow-band filter 
the combined narrow-band image is about 1 magnitude shal- 
lower than the combined U-band image for objects with flat 
SEDs. We can therefore only detect a deficit of emission in 
the narrow filter for relatively bright U-band selected objects 
(brighter than U(AB)=25.3). As candidates we select galaxies 
with colours above n(AB)— U(AB)=0.6 within la. We detect 
two and five such galaxies in the field of GRB 00030 1C and 
GRB 000926 respectively. Their positions are indicated with 
squares in Fig. ||. Deep follow-up spectroscopy is required to 
confirm if these objects are high redshift galaxies with strong 
z=2.04 intrinsic or intervening Lya absorption lines. 

3.6. The host galaxies 

We now discuss what can be inferred on the properties of the 
host galaxies from the observations presented here. The host 
galaxy of GRB 000301C is not detected in any of the bands. 
In Fig. H we show a 10x10 arcsec 2 region centred on the po- 
sition of the afterglow of GRB 00030 1C from the combined I- 
band image and from a weighted sum of the combined Lya and 
U-band images. We calculate the weight as the inverse of the 
variance of the background after scaling to the same flux level. 
In the I-band we detect the red galaxy 2 arcsec north of the af- 
terglow also seen in the STIS images by Smette et al. (2001). 
In the Lya + U-band image we detect no significant emission 
from either this red galaxy or from the host of GRB 00030 1C 
found by Fruchter et al. (2001a) above an estimated 2a limit 
of U«27.7 per arcsec 2 . Since the host is very faint (R=28-30, 
Bloom et al. 2002; Levan private communication; Fruchter et 
al. in prep) this non-detection is not surprising. 

The host galaxy of GRB 000926 is detected in all bands 
and is a strong Lya emitter with a rest-frame equivalent width 
of the Lya emission line of 7 1 A. The galaxy is shown in more 
detail in Fig. ^. Deep HST images of the host has been reported 
by Price et al. (2001) and Castro et al. (2001). These images 
confirm what was found from the ground based R-band images 
(Fynbo et al. 2001b and Fig. ^), namely that there are several 
compact objects within a few arcsec, one of which coincide 
with the position of the OT. To better compare the extent of the 
Lya emission to that of the continuum we retrieved and driz- 
zled (Fruchter & Hook 2002) the May 19 2001 F606W HST 
imaging data of the host from the HST archive. These data are 
published in Castro et al. (2001). In the top panel of Fig. ^ we 
show a section of the F606W image together with the contours 
of the Lya emission. The Lya emission covers all of the com- 
ponents seen in the continuum and we conclude that they are 
all at the same redshift and hence part of the host galaxy and/or 
its environment. As a simple model we fit the Lya image with 
two point sources using DAOPHOT (Stetson 1987, 1997). The 
derived distance between the two best fitting point sources is 
2.22±0.08 arcsec (corresponding to 20.0 kpc), which is consis- 
tent with the distance between the brightest knots in the HST 
image (2.25±0.01 arcsec). From the point source fit we find 
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GRB000926: Lya GRB000926: R-band 
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GRB000926: U-band GRB000926: l-band 




Fig. 6. Upper plot: The position of GRB 00030 1C as imaged 
in the I-band (right) and Lya + U-band (left). The size of the 
images is 10x10 arcsec 2 . No significant emission at the posi- 
tion of the afterglow is detected to a 2<r limit of U(AB)=27.7 
per arcsec 2 . Lower plot: The host galaxy of GRB 000926 as 
imaged in Lya, U, R and I. The R-band image is taken from 
Fynbo et al. (2001b). East is to the left and north is up. The 
size of the images is 10x10 arcsec 2 . The position of the opti- 
cal afterglow is indicated with an x . 



that about 65% of the Lya emission is emitted from the east- 
ern compact knot and 35% from the western knot. The position 
of the western (Lya) knot is spatially coincident with the po- 
sition of the OT to within 2a of our error of about 0.10 arcsec 
on the relative astrometry between the R-band image with the 
OT from Fynbo et al. (2001a) and the Lya image. Castro et al. 
(2001) measure an impact parameter of the OT relative to the 
centroid of the western (continuum) knot of 32±4 mas based 
on five epochs of F606W images of the OT superimposed on 
the host galaxy. 

As seen from Fig. ||the western and eastern knots seem to 
have different broad band colours with the western component 
being the reddest. To better quantify this we also retrieved the 
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Table 2. Photometric properties of the 12 candidate LEGOs of rank 1 in the fields of GRB 000301C and GRB 000926. Aperture 
corrections for the 3'.'0(4'.'0) aperture are 0.36(0.21), 0.27(0.15) and 0.23(0.14) mag for narrow, U and I respectively. 



Object 


n(AB) 


U(AB) 


I(AB) 


f(Lya) 

10" 17 erg s" 1 cm" 2 


W cm (obs) 


SFR-Lya 












A 


M yr" 1 


Isophote 
















S301.1 


25.07±o?7 


25.981°;- 


24 


70+0.13 

• Z8 -0.11 


3.36± 0.54 


69+H 


0.92±0.15 


S301.2 


24.88l°i* 


26.161°;?! 


26 


Q ,1 + 1.66 

•^-0.63 


3.98± 0.58 


13011? 


1.09±0.16 


S301.3 


24.98+°;i« 


26.80l°;It 


26 


^7+1.89 

■ JA -0.65 


3.64± 0.52 


- , - £ -- J _144 


1.00±0.14 


S301.4 a 


23.9llg-.io 


25.1312;!? 


25 


44+0-75 


9.80± 0.87 


120+ 3 ? 


2.68±0.24 


S301.5 


24.54l°;i 3 


25.501°;" 


24 


fin +0.25 
SU -0.20 


5.46± 0.64 


7 c+28 
/; >-22 


1.49±0.18 


S301.6 11 


25.481°;- 


26.90l°;« 


25 


17+0.24 
- 1 '-0.19 


2.30± 0.42 


163+ 7 f 


0.63±0.11 


S301.7 


70 1 i +0.06 
AJ - L i -0.06 


24.87lg;io 


25 


67+1.41 

■ D '-0.59 


20.40± 1.05 


292+g 


5.58±0.29 


S301.8 


24 n+°- n 


25.18l°;i 


23 


77+O.IO 
■ /z -0.09 


7.95± 0.75 


881?? 


2.18±0.21 


S301.9 


23.30±°:°S 


94 Q0+ 011 


24 


«+0.31 
■ JO -0.24 


17.14± 1.00 


223l|5 


4.69±0.27 


i'.'O aperture 
















S301.1 


24.58l°; 2 2 ? 


95 4.Q+ - 22 


23 


oo+0. 17 
• S8 -0.15 


5.27± 1.04 


- 


1.44±0.28 


S301.2 


24.29+°,-?? 


9c 77+O.I8 
zj.z/_ 15 


25 


Q7+3.90 
■ yA -0.7i 


6.89± 1.06 


- 


1.88±0.29 


S301.3 


24.32+°-'?? 


26 15+ - 45 

ZO.lJ_ 32 

95 07 +0 - 15 
zo.u/_ 13 

2 C 1A+0.15 


25 


O7+3.90 
■ yA -0.7i 


6.72± 1.06 


- 


1.84±0.29 


S 301.4° 


+0.11 
ZJ./f_ 10 


25 


07+2. 36 
■ Oz -0.69 


11.39± 1.09 


- 


3.12±0.30 


S301.5 


24 16+°- 17 

Zt.lO_0.15 


24 


or,+0.27 
JKJ -0.21 


7.73± 1.07 


- 


2.12±0.29 


S301_6 a 


24 41+0-23 


9C C9+0.23 
iJJZ -0.19 


23 


QA+0.19 

■ yH -0.16 


6.05± 1.05 


- 


1.66±0.29 


S301.7 


70 i i +0.06 
"• x -0.06 


24.83lg;i§ 


25 


O7+3.90 
■ yA -0.7i 


20.47± 1.11 


- 


5.60±0.30 


S301.8 


23.87l°; 1 12 : 


24.86l°i? 


23 


on+0.11 
■ Jy -0.10 


10.14± 1.08 


- 


2.77±0.30 


S301.9 


11 97+0.07 
zj.z/_ 07 


24.90lg;i? 


24 


64+0-39 


17.54± 1.10 


- 


4.80±0.30 


4!0 aperture 
















S301.1 


24.59l°|* 


25 27+ 24 


23 


ni+0.24 
■ yL -0.20 


5.20± 1.26 


- 


1.42±0.34 


S301.2 


24.331 ,; 2 ? 


25.09+°-?? 


25 


69+ 4 -20 

■ uz -0.74 


6.64± 1.28 


- 


1.82±0.35 


S301.3 


24.141°^ 


25.94l°jl 


25 


62 +4 ' 20 

■° -0.74 


7.87± 1.30 


- 


2.16±0.36 


S301.4 a 


23.591°;!? 


24 73+ - 14 


25 


69+ 4 - 20 
■ DZ -0.74 


13.10± 1.33 


- 


3.59±0.37 


S301.5 


24.081°;?? 


95 14+0- 21 


24 


m+0.27 
ui -0.22 


8.33± 1.30 


- 


2.28±0.36 


S301-6 a 


94 94+0-23 
Z ^- ZH -0.19 


25. 10+°-?? 


23 


IO+0.11 

lu -0.10 


7.23± 1.29 


- 


1.98±0.35 


S301.7 


23 02+°' 07 


24.771°;^ 


25 


' u:, -0.74 


22.09± 1.36 


- 


6.05±0.37 


S301.8 


2 t 75+0.14 

z ->-' J -0.13 


24.771°;^ 


23 


,17+0.15 
' — 0.13 


11.26± 1.33 


- 


3.08±0.36 


S301.9 


91 T 5+0-08 
ZJ>.1J_ Q7 


24.841°,;!° 


24 


O/1+0.70 
••"--0.42 


19.75± 1.36 


- 


5.41±0.37 


Isophote 
















S926_l ft 


91 15+0-07 
"■JJ_0 .07 


94 Q4.+ 013 


24 


17+0.15 
1 '-0.13 


16.30± 1.02 


9 1 7+62 
zl/ -45 


4.46±0.28 


S926.2 


22.78l°-°4 


+0.10 
Zt.JZ_ Q g 


25 


oc+1.21 
■ZJ-0.56 


27.73± 1.14 


287l° 2 


7.59±0.31 


S926.3 a 


25.501°;?! 


27 04+ 042 
Z/.Ut_ 30 


24 


7Q+0.10 
•'-'-0.09 


2.27± 0.40 


204! 2 ! 8 


0.62±0.11 


3f. aperture 
















S926.1 6 


23.48l°;°« 


94 Q4+ 013 

Zt.5"t_ 12 


24 


7 4 +0.17 
■ ZH -0.14 


14.49± 1.04 




3.97±0.28 


S926.2 


22.83l°-gt 


24 60+°- 10 

zt.OU_ g 


25 


C7+0.72 
■-"-0.43 


26.39± 1.06 




7.22±0.29 


S926.3 a 


24.9ii°;i? 


25.66l°; 22 


22 


, 6 +0.02 
io -0.02 


3.88± 0.97 




1.06±0.26 


4'.'0 aperture 
















S926.1 6 


91 91+0.09 
ZJ.ZJ_ 08 


24.751°;^ 


23 


O1+0.15 
■ 81 -0.13 


18.26± 1.38 




5.00±0.38 


S926.2 


22.701°;°! 


24.57lg;i? 


25 


oq+2.10 
Oy -0.67 


29.67± 1.40 




8.12±0.38 


S926.3 a 


24.68l°; 3 2 « 


25.60l°; 3 2 ? 


21 


60 +0.02 
■ DU -0.02 


4.80± 1.27 




1.31 ±0.35 



a Broad-band photometry affected by bright nearby object. 
b The host galaxy of GRB 000926. 



F814W HST image from May 20 (lower panel in Fig. [7]- not 
drizzled due to the lower S/N of the host in the F814W im- 
age) and measured the magnitudes of the two brightest compo- 
nents of the host in F606W and F8 14W images. In these images 
we measure the AB colours (V606 — Ism) as = —0.11 ±0.11 
and (V606 — 1&ia)ab = 0.19±0.07 for the eastern and west- 
ern components respectively. Using the relation between (V606 
— ^%ia)ab, redshift and the spectral slope (3 (fx oc A^) de- 
rived by Meurer et al. (1999) we derive spectral slopes (3 = 
— 2.4±0.3 and — 1.4±0.2 for the two components. Based on 



the relation Aiqqq = A A3 ± 1.99/3 also derived by Meurer et 
al. (1999) we conclude that the eastern component must be es- 
sentially dust free (the derived value of A1600 is formally nega- 
tive, but consistent with 0), whereas for the western component 
we find ^4i6oo = 1.6±0.4. The luminosity-weighted mean value 
of A1600 for the sample of Lyman-Break galaxies studied by 
Meurer et al. (1999) is 1.8±0.2. 

The restframe EW of 711?° A is within the range of 50- 
200 A expected for dust-free galaxies (Chariot & Fall 1993). 
Furthermore, the SFR derived from the Lya flux in the 4'.'0 
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Table 3. Photometric properties of the 7 candidate LEGOs of rank 2 in the fields of GRB 000301C and GRB 000926. Aperture 
corrections for the 3'/0(4'/0) aperture are 0.36(0.21), 0.27(0.15) and 0.23(0.14) mag for narrow, U and I respectively. 



Object 


n(AB) 


U(AB) 


I(AB) 


f(Lya) 
10" 17 erg s _1 cm -2 


W cm (obs) 

A 


SFR Lya 
Mq yr" 1 


Isophote 














S301.10 a 


25.4018;?! 


26.57l8;l 4 


21 67 +0 ' 01 


2.48± 0.42 


I08ljj 


0.68±0.11 


S301.ll 


25.7418;1I 


26.6218;?! 


27.0218;?? 


1.81± 0.36 


65l 4 3 l 


0.49±0.10 


y.'O aperture 














S301.10 a 


24.67lg : g 


9c co+0.24 
iJ ' JO -0.20 


2o.54i8.81 


4.83± 1.03 


- 


1.32±0.28 


S301.ll 


24.55t° : * 


25.4918;?! 


25.0718;^ 


5.41± 1.04 


- 


1.48±0.29 


4'.'0 aperture 














S301.10" 


24.6St° fs 


25 74+ - 40 
zj./t_ 2g 


20.4318:81 


4.79± 1.25 


- 


1.31±0.34 


S301.ll 


94 4.O+0-30 
"■*"-0 .24 


9S 94+0-23 


24 7Q+0-66 


5.72± 1.26 


- 


1.57±0.35 


Isophote 














S926.4 


9S n4 +0 - 21 


25.7618;!! 


24.0618:8? 


3.46± 0.56 


47111 


0.95±0.15 


S926.5 


9^ tq+0.22 
■iJ-JJ-0 .18 


26.2518;?^ 


24.78l8i? 


2.64± 0.44 


70l 4 3 4 o 


0.72±0.12 


S926.6 


24.72l°;i« 


25.4518:11 


9Q 71 +0.06 
z.j. ; i_o.o6 


4.64± 0.65 


48l? 4 


1.27±0.18 


S926.7 


24.7918 It 




23.8918.8^ 


4.34± 0.57 


52l? 4 


1.19±0.16 


S926.8 


25.64lgiH 


26.5918;!? 


25.8oi8;18 


1.98± 0.39 


74lf 7 


0.54±0.11 


3'.Q aperture 














S926.4 


24.4618;?! 


25.nl8;| 


23.74l8;J8 


5.91± 1.00 




1.62±0.27 


S926.5 


24.8ll8;g 


25.65l 22 


23.94l8;l? 


4.28± 0.97 




1.17±0.27 


S926.6 


94 97+0-18 


25 14+ - 16 


29 4fr+ 003 


7.00± 1.01 




1.91±0.28 


S926.7 


24.48t°;S 


24.911811 


23.4818.8? 


5.76± 0.99 




1.58±0.27 


S926.8 


24?t -+0.30 


25.6ll8;H 


24.95l8;i 


4.50± 0.98 




1.23±0.27 


4'.'0 aperture 














S926.4 


24.39lg;ll 


24.97l8;18 


23.62+8'!? 


6.25± 1.30 




1.71±0.36 


S926.5 


24.83l8; 4 3 l 


9S SQ+0-36 

z.j.J7_q 27 


23.951°^ 


4.19± 1.26 




1.15±0.34 


S926.6 


24.2ll°;S 


94 01 + 18 


91 04+0.02 
■"■"-o .02 


7.37± 1.32 




2.02±0.36 


S926.7 


24.4818:11 


24 R4+ - 17 


90 9f\+0.08 

Z.J.Z.W_g gg 


5.75± 1.29 




1.57±0.35 


S926.8 


24.60l8;| 


25.881°;" 


94 £53+0.30 


5.20± 1.28 




1.42±0.35 



a Broad-band photometry affected by bright nearby object. 



aperture taking into account the aperture correction is 6 Mq 
yr~ 1 . This is only about a factor of 2 lower than what we found 
based on the restframe 2100 A continuum flux (Fynbo et al. 
2001a). This means that there cannot be large amounts of dust 
in the host galaxy to absorb the resonantly scattered Lya pho- 
tons. This is consistent with the low values of A1600 derived 
above and with the low SMC-like extinction measured from 
the optical afterglow (Fynbo et al. 2001b; Price et al. 2001). 

Of all the candidate LEGOs, the Lya emission from the 
host galaxy of GRB 000926 is the only to have an obvious 
multi-component morphology. This suggests that we see sev- 
eral galaxy sub-clumps (or LEGOs as we call them) in the pro- 
cess of merging and that the GRB could be related to a merging 
induced starburst. 

If we assume a R-band magnitude of R=28 for the 
GRB 00030 1C host galaxy it is about 40 times fainter than 
GRB 000926 host galaxy. If the GRB 00030 1C host has the 
same colour and Lya equivalent width, its Lya flux will be 
^4x 10~ 18 erg s _1 cm~ 2 which is a factor of 5 below the (5(t) 
detection limit in the narrow-band image. Hence, our data does 
not exclude that the GRB 00030 1C host galaxy is also a high 
equivalent width Lya emitter. 



3. 7. Exclusion of contamination due to foreground 
emission line sources 

For searches for high redshift emission-line galaxies it is im- 
portant to realize that lower redshift galaxies with other emis- 
sion lines in the narrow filter can mimic high redshift Lya can- 
didates (e.g. Stern et al. 2000; Fynbo et al. 2001c). For the host 
galaxy of GRB 000926 the redshift is known from the spec- 
troscopy of the Optical Transient (OT) (Fynbo et al. 2001d; 
Castro et al. 2001; M0ller et al. in prep.). For the remaining 
candidate LEGOs we here consider other possible emission 
line sources than Lya at z=2.04. We can exclude foreground 
galaxies with [Oil] in the narrow filter, since such objects 
would have to be at redshift z=0 and hence be much brighter 
than our candidates. Distant planetary nebulae in the galaxy 
could be selected due to their strong O II emission lines, but 
such planetary nebulae are very rare at the high galactic lat- 
itudes of our fields (+44.4° and +37.3° respectively). In total 
only about 10 planetary nebulae are known in the Milky Way 
halo (Howard et al. 1997). Other possible explanations are C IV 
1550 A at z= 1.39 or Mg II 2800 A at z=0.32. The presence of 
strong C IV and Mg II emission requires an underlying AGN. 
The surface density of AGN at z<2.2 down to a B-band mag- 
nitude of 22.5 is 130 deg~ 2 (Hartwick & Schade 1990). If we 
extrapolate to B=25.0 using the observed faint end slope of 
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Fig. 7. Top panel: The WFPC2/F606W drizzled image of the 
host galaxy of GRB 000926 from Castro et al. (2001) plotted 
together with the contours (shown with dotted lines) of the Lya 
emission. North is up and east is to the left. The scale of the 
image is indicated with a bar in the upper left corner and the 
estimated error on the alignment of the HST and NOT images 
is shown in the lower left corner. The Lya emission covers all 
the components of the continuum emission showing that they 
are all part of the host galaxy and/or its environment. Most 
of the Lya emission (about 65%) is emitted from the eastern 
most compact knot seen in the F606W image. Lower panel: 
The same region as shown above, but from the (undrizzled) 
F814W image. 



d\ogN(B)/dB = 0.35 we find 800 deg" 2 . Hence, we expect 
6 z<2.2 AGNs to this limit in each of our 27.6 arcmin 2 fields. 
The probability that these fall in the small volumes at z=0.32 
and z=1.39 (0.18% and 1.9% of the total comoving volume out 
to z=2.2 respectively) probed by our filter is small, especially 
since the space density of AGNs is declining rapidly with de- 
creasing redshift. We therefore do not expect a significant C IV 
and Mg II contamination due to AGN. In conclusion, the only 
probable explanation for the excess emission in the narrow fil- 
ter is Lya at z=2.04. 



4. Discussion and conclusions 

4.1. LEGOs 

We first compare the photometric and other properties of 
the candidate LEGOs in the fields of GRB 00030 1C and 
GRB 000926 to properties of z«2 galaxies found in other sur- 
veys. 

Fontana et al. (2000) have determined photometric redshifts 
for galaxies in the HDF and NTT deep fields. From their cat- 
alogue we find 22 galaxies with photometric redshifts within 
Az=0.15 of z=2.04. The I(AB) magnitudes of these galaxies 
range from I(AB)=23.7 to I(AB)=25.8. The range of I(AB) 
magnitudes for the candidate LEGOs in our fields is extend- 
ing to somewhat fainter magnitudes as 4 of the 19 candidates 
are not detected above the 2a level of 25.6 in the 3'.'0 aperture. 

Kurk et al. (2000) and Pentericci et al. (2000) have de- 
tected and spectroscopically confirmed 14 LEGOs in the field 
around the radio galaxy PKS1 138-262 using narrow-band ob- 
servations. They reach a 5<r detection limit of 1.4 x 10~ 17 erg 
s _1 cm~ 2 over a field of 38 arcmin 2 at z=2.156. Their filter 
spans a redshift range of Az = 0.054. Their survey parameters 
are therefore similar to ours. The inferred number of LEGOs 
per arcmin 2 per unit redshift in their survey is 6.8±1.8. We ex- 
pect to confirm 100% (75%) of our rank 1 (rank 2) candidates, 
hence 17 LEGOs total in two 27.6 arcmin 2 fields observed with 
a filter spanning Az = 0.037. This corresponds to 8.3±2.0, 
which is consistent with the density found by Pentericci et al. 
(2000). 

Pentericci et al. (2000) estimate that their field is a factor of 
6 overdense relative to the field, but this estimate is very uncer- 
tain as it is based on ;) only the four brightest objects in their 
sample, ii) a comparison with the spike of Lyman-Break galax- 
ies observed by Steidel et al. (2000) at a significantly larger 
redshift (z=3.09) and with different selection criteria, and in) an 
assumption about the (unknown) fraction of galaxies at z=2. 16 
that are Lya emitters. Hence, we consider the evidence for an 
overdensity to be at best weak. Unfortunately, there are cur- 
rently no similar (in redshift and depth) surveys for Lya emit- 
ters in blank fields so we do not know yet whether the QSO 
and GRB fields are both equally overdense or if they represent 
the mean density of LEGOs at z«2 (which we consider more 
likely). We are currently undertaking a large blank field sur- 
vey at the Nordic Optical Telescope using the same filter with 
the purpose of examining the density and spatial distribution of 
LEGOs in blank fields. 

4.2. GRB host galaxies 

This study confirms what is already known from broad-band 
studies, namely that GRBs occur in galaxies covering a very 
broad range of luminosities. The faintest host galaxies de- 
tected so far in the broad-bands are the hosts of GRB 00030 1C 
and GRB 990510 at R«29.7 and R«29.2 and the host 
of GRB 980326 at V?»30.3 (Levan, private communication; 
Fruchter et al. in prep.). 

We have used deep Lya narrow band imaging to map out 
the star formation in volumes around the host galaxies of two 
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GRBs. GRB 000926 occurred in one of the strongest centres 
of star formation within several Mpc, whereas GRB 00030 1C 
occurred in an intrinsically very faint galaxy far from being the 
strongest centre of star formation in its galactic environment 
(cf. Table g). 

What is the underlying mechanism for GRBs? The answer 
to this question is still not known, but we now have several vital 
clues. It has been noted that GRB host galaxies generally are 
UV bright and have emission lines with large equivalent widths 
(Bloom et al. 2002; Fruchter et al. in prep.), indicating that they 
are actively star forming. This result is supported by our detec- 
tion of strong Lya emission from the host of GRB 000926, 
and is not contradicted by the lack of Lya emission from the 
host of GRB 00030 1C where even a very high Lya rest equiv- 
alent width of ~450 A would have remained undetectable in 
our deep narrow band image. It is quite natural therefore, to 
hypothesize that the GRB phenomenon is linked to star forma- 
tion (Paczyhski 1998). A strong prediction from this hypothe- 
sis is that the GRB host galaxies are drawn from the underly- 
ing population of high redshift galaxies weighted with the in- 
tegrated star formation rate in each magnitude bin. If the LF of 
the underlying population was flat, then we would predict that 
all GRB host galaxies should be very bright. Conversely a very 
steep LF would result in all GRB host galaxies being drawn 
from the very faintest galaxies. The observation that GRB host 
galaxies span a very wide range of luminosities therefore places 
constraints on the allowed range of faint end slopes of the LF. 
For example, if we use a standard Schechter form 



<$>(L)dL = $* 



exp 



L 



of the LF, and use R*=23.50 then a simple calculation shows 
that the observed flat distribution in the range R=24 to R=30 
requires a faint end slope a = —1.9 ±0.7. Such a steep 
(steeper than in the local universe) faint end slope is in good 
agreement with the faint end slope for Lyman-Break galaxies 
at z=3.0 (a = -1.57 ± 0.11, Adelberger & Steidel, 2000), and 
with the faintness and small impact parameters of high redshift 
DLA galaxies (Fynbo et al. 1999; M0ller et al. 2001). The cur- 
rent data set is therefore compatible with the hypothesis that 
GRBs trace star formation, and the outlook for impending ad- 
vance on this question is good. A full statistical treatment of a 
slightly larger sample of GRB host galaxies would significantly 
decrease the error bars on a, and would show if it remains com- 
patible with the faint end slope of the LBGs. More importantly, 
the GRB host galaxies probe fainter magnitudes than the LBG 
samples and will therefore reveal if there is a turn-over, or if the 
faint end slope continues. The question of how far the faint end 
of the LF continues with unchanged slope is important for the 
determination of the total star formation rates at high redshifts. 
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Appendix A: Errors on colours for faint sources 

Here we detail a method by which to calculate two-sided error- 
bars for sources detected at low signal-to-noise ratio. We wish 
to calculate the colour m(a) — m(b) for a source with fluxes a 
and b in two passbands. We first derive the likelihood function 
for the flux ratio a = b/a. To do this we introduce (3 = a ■ b 
and note that b 2 = a ■ (3 and a 2 — (3 /a. The joint likelihood 
function for a and (3 can then be found as 



L(a, f3) oc cxp (- 



where (cto, /3rj) is the maximum likelihood estimate. By maxi- 
mizing the likelihood function in the variable (3 for a fixed value 
of a we find 



bo 



an 
a 'i a 



1 



and 



1 

aer- 



After inserting these expressions we find the likelihood func- 
tion for a of the form 



logL(a) = logL(ao) 



1 ala 2 a {a - a ) 2 + &§o£(l/a - l/a ) 2 



We then determine the two-sided error-bars for the ratio a as 
the two solutions to the equation 

logL(a + a a) = logL(ao) - 1/2 

Finally, we use the error-bars on a to calculate the error-bars on 
the colour in the same way as for the magnitudes in Sect. 3.2. 
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Fig. 5. Regions of size 10 x 10 arcsec 2 around each of the eleven candidates LEGOs in the field of GRB 00030 1C and the eight 
candidates LEGOs in the field of GRB 000926. For each candidate we show a sub-image from all three filters narrow (top row), 
U (middle row) and I (bottom row). The name and rank of the candidates are printed above each of the narrow-band sub-images. 
The host galaxy of GRB 000926 is S926_l . 



